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Abstract 

Recent a s determinations from jet production at HERA are presented. Three 
different observables lead to consistent a s values. The most accurate result is: 
a s (M z ) = 0.1181 ±0.0030 (exp.)±°;°°^ (thco.)i[™ (pdf) 
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1. Introduction 

Perturbative Quantum Chromodynamics (pQCD), 
the theory describing short distance strong 
interactions, has a single free parameter, the strong 
coupling constant a s . This coupling depends on 
the renormalisation scheme and the energy scale. 
At the reference scale, which is customarily chosen 
as the mass of the Z° boson, the strong coupling 
is presently known to an accuracy of about 4% JD . 
This is more than a factor of 10 less precise than 
the coupling of the weak force. Measurements of 
a s over a wide energy range provide a fundamental 
test of the evolution predicted by the QCD (3- 
function and can severely constrain physics beyond 
the Standard Model. 

Precise determinations of a s have been per- 
formed in e + e~ annihilation processes over a large 
range of centre of mass energies from 14 < yfs < 189 
GeV using jet rates, event shapes, r decays and the 
inclusive ratio of hadron to lepton branchings H . 

Jet production at hadron colliders is a 
particularly well suited process to measure the 
energy dependence of a s in one single experiment. 
However, the presence of hadrons in the initial state 
leads to the complication that the measured cross 
section always depends on the product of a s and 
the parton density functions (PDF). 

The high centre of mass energy (s) of HERA, 
colliding 27.5 GeV positrons on 820 GeV protons, 
allows hard processes to be studied in deep-inelastic 
scattering (DIS). Particularly well suited for pQCD 
studies is the Breit frame in which the photon 
is purely space-like and collides head-on with the 
proton. Besides the photon virtuality, Q 2 , the 
transverse energy Et produced in the Breit frame 
provides a hard scale. 

Jet cross sections at high Et are collinear and 
infrared safe observables which can be calculated in 
pQCD. 
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Figure 1. Examples of 0(a s ) and O(o%) Feynman 
diagrams dijet events in e p -collisions. 

The jet cross section cr jet can be expanded to: 
y] «?(/*?) / -T C a , n (x Bj /£, IJ, 2 r , n), ...) f a/p {f, Ll 2 f ), 

where x Bj = Q 2 /{2p ■ 7), £ = x Bj (1 + Mf^/Q 2 ) 
is in leading order (LO) the fractional momentum 
of the struck parton in the proton and Mjj the 
invariant dijet mass. \i r (/x/) is the renormalisation 
(factorisation) scale. The perturbatively calculable 
coefficients C a , n are presently known to next-to- 
leading order (NLO). The jet cross section is 
directly sensitive to a s and the gluon (f g / p ) and 
quark (f q / p ) densities in the proton. For low Q 2 
(below rs 40 GeV 2 ) jet cross sections calculated in 
NLO strongly depend on fi r [B|||. pQCD analyses 
are therefore better performed at high Q 2 , where 
also hadronisation effects are smaller H . 

Double differential inclusive dijet cross sections 
d 2 a/d^dQ 2 and d 2 a/dxdQ 2 for 200 < Q 2 < 5000 
GeV 2 , where jets are defined by the inclusive Kt 
algorithm j| asking for E T > 5 GeV , E T ,i + 
Et.2 > 17 GeV and -1 < r)i a b < 2.5 have 
been used for a preliminary extraction of the gluon 
density in the range 0.01 < £ < 0.1 at fixed ^ — 

200 GeV 2 ||. In this analysis a s has been set to 

f The pseudo-rapidity is defined by r\ = In tan 9/2. 



the world average. F 2 data || for 200 < Q 2 < 650 
GeV 2 have been included in a combined fit to get 
an additional constraint on the quark densities. 

At this conference new preliminary results have 
been presented where the opposite strategy has 
been pursued. The parton densities have been taken 
from global fits to inclusive DfS data and exclusive 
high--Er processes and then a s is fitted. Either sets 
of PDF with different input a s have been used to 
account for the anti-correlation between a s and the 
PDF or the dependence on the PDF is included in 
the systematic error. The final aim is to perform a 
consistent simultaneous fit of both a s and the PDF. 

2. a s from Inclusive Jet Cross Sections 

The single inclusive jet cross section da/dQ 2 dET 
has been measured in four different Q 2 regions 
within 150 < Q 2 < 5000 GeV 2 and 0.2 < y = 
Q 2 /(s Xbj) < 0.6 by the HI collaboration. Jets 
are defined by the inclusive Kt algorithm requiring 
E T > 5 GeV in the range -1 < r\ lab < 2.5. The 
data are corrected for detector effects. Over the 
whole phase space reaching Et up to 50 GeV a 
NLO calculation based on DfSENT H] using the 
CTEQ5M § PDF with /z 2 . = E T and ^ 2 = 200 

GeV 2 describes the data. Hadronisation effects 
lower the NLO prediction by about 3-10% §. For 
each Et and Q 2 bin a s is adjusted to the measured 
cross section using a \ 2 minimisation where the 
systematic errors are included in the fit M. All 
results are consistent with each other. In Fig. @ the 
results are shown as a function of Et after having 
combined the four Q 2 bins. The fitted a s evolve as 
predicted by the renormalisation group equation. 
The combined a s (M§) result is: 
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0.1181 ± 0.0030 (exp.) +0;°°*j (theo.) + 00036 



-0.0015 



(pdf.) 



The largest contribution to the experimental error 
comes from the uncertainty on the hadronic energy 
scale. The theoretical error is mainly determined 
by the uncertainty of the hadronisation corrections 
and by the \x r dependence. The contribution 
obtained by varying /i 2 by a factor of 4 amounts 
to io'oo34- The hadronisation uncertainties are 
estimated in a very conservative way using the QCD 
Monte Carlo models. Since the spread between 
this models is rather small, half of the size of 
the hadronisation corrections, but at least 3%, is 
assigned as systematic uncertainty. If ^jl 2 = Q 2 
instead of [i 2 — E T is used a consistent a s {M^) 
result is obtained: 

0.1221 ± 0.0034 (exp.) +°;°™$ (theo.) t°;°oi1 (pdf) 

For /i 2 = Q 2 the fx r dependence increases by 
« 40%. 
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Figure 2. a s fit result versus the jet Et- The lines 
indicate the prediction of the renormalisation group 
equation for the combined results and their uncertainty. 

The systematic uncertainty on the PDF is 
difficult to estimate, because none of the global 
analyses has yet provided PDF errors. Since 
the fit results of the different groups, MRST @ 
and CTEQ [||, are based on similar data 
sets and assumptions, the spread of different 
parametrisation does not give a realistic error 
evaluation. Recently the CTEQ collaboration [O 
has investigated possible variations on the gluon 
distribution using different parametrisations and 
allowing the quality of the fit to be degraded. The 
resulting parametrisation is the first step towards 
the estimation of uncertainties, but can not replace 
a proper error analysis. 

Another point of concern is the influence of 
the a s value used when deriving the PDF. If a 
strong correlation between the fit result and the 
initial assumption on a s was found, it would be 
questionable whether the ansatz to fix the PDF 
and to extract a s would be meaningful. This 
dependence can be tested by means of series of 
parton distributions scanning from low to high 
values of a s . 

A comprehensive study of the fitted a s on all 
available recent PDF, is presented in Fig. 0. The 
largest deviation to the central result obtained 
with the CTEQ5M parametrisation is quoted as 
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input a s (M z ) used in global fit of parton densities 
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Figure 3. Dependence of the a a fit results on the 
assumed PDF in four Et regions. 

systematic error. Only for the highest Et a slight 
dependence on the a s (M|) value used in the global 
fits is found. For all other bins no correlation is 
observed. This is different from the observation 
made by the CDF collaboration pjj in an analysis 
of single inclusive jet cross sections measured in pp 
collisions where a strong correlation between the a s 
used in the PDF and the fit result was found. 



3. a s from Dijet Rates 

The ZEUS collaboration has measured dijet cross 
sections and rates in the phase space 470 < Q 2 < 
20000 GeV 2 and < y < 1. Jets are defined as 
in the dijet analysis described before, except that 
the two highest Et (and not all) jets have to fulfil 
— 1 < r/iab < 2. In addition, only events with 
exactly two jets are considered. The exclusive dijet 
cross sections as a function of Et, Mjj, r\i a b, xbj , 
Zp.i = Ei {1-cosQi) I \Y,k=x,z E k( l - cos °k)), and£ 
are - both in shape and in absolute magnitude - well 
described by a NLO calculation using CTEQ4M p| 
for /i 2 = fj,% = Q 2 . The data are corrected for 
detector, electroweak radiation and hadronisation 
effects. Hadronisations corrections are at most 
10% and decrease with increasing Q 2 . For Q 2 > 
5000 GeV 2 the contributions from Z° exchange 
become visible. They have been subtracted from 
the data. 

The good description of the data allows an ex- 
traction of a s fromi? 2 +i(Q 2 ) = V2+i(Q 2 ) / o- to t(Q 2 ) ■ 
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Figure 4. Exclusive dijet rate as a function of Q 2 . 

The dijet rates shown in Fig. [f| are compared to 
NLO QCD calculations based on the CTEQ4A PDF 
series which have been fitted for different a s as- 
sumptions. This series is also used to parametrise 
the anti-correlation between a s and the PDF G3. 
The value of a s is obtained from a % 2 minimisation: 



0.120 ±0.003 (stat. 



■ ±8:885 («p0: 



:8:88i (th«.) 



The experimental systematic error is dominated 
by the uncertainty on the hadronic energy scale 
(±0.005). QCD Model dependencies to unfold the 
data to parton level account for an error of io^- 
The theoretical uncertainty is estimated by varying 
the renormalisation scale /j, 2 — Q 2 by a factor of 2 
(±0.002), by replacing the assumed PDF by MRST 
sets (-0.002) and by taking DISASTER++ [HI 



instead of DISENT (±0.002). If (if is varied by 
a factor of 4 as in the previous analysis, the /i r 
uncertainty increases to ±0.004. The dijet rates 
change only by at most 2.5% when the PDF from 
the CTEQ gluon uncertainty study is used. The 
influence on the a s result is therefore expected to 
be small. 

4. a s from Differential Jet Rates 

HI has measured l/N to t {dn/dQ 2 dy2) for Q 2 > 
150 GeV 2 and 0.1 < y < 0.7. All particles are 
clustered on the basis of a distance measure y 
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j j /scale, where di t3 



2mm(Ef,E 2 )(l- cos 9ij), 
1 jets are found . j/2 is the smallest value 



f +1 denotes the proton remnant. Ei is the energy of the 
particles and 8ij the angle between them. 



of y where two jets can be resolved. It is therefore 
a measure for the jet structure. Two jet algorithms 
are used. In both cases the data are described by a 
NLO calculation based on CTEQ5MJ0. 

The Kt jet algorithm for DIS [|16| is applied 
in the Breit frame and scale is set to 100 GeV 2 . 
The proton remnant p is included as particle with 
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2E 2 (1 - cos6i). After a cut y 2 > 0. 



only events with Et > 10 GeV in the Breit frame 
remain. A x 2 minimisation gives a s (Af§): 



n 1 143 +00075 / s +0.0074 /,, s +0.0008 

U.J.110 _ .0089 l LX P-J -0.0064 l+UCO.j _ .0054 



(pdf) 



A variant of the JADE algorithm is applied 
in the HERA laboratory frame and the squared 
invariant mass of the hadronic final state W 2 is 
used as scale, applied in the HERA laboratory 
("mod. Durham algorithm"). A massless particle 
is added in the clustering procedure to account for 
the escaped longitudinal momentum. y 2 is required 
to be above 0.005 to select hard processes. The fit 
result is: 

fl 1 1 SQ +0.0064 / \ +0.0059 (+y, Dr , A +0.0013 /_ jr \ 

U.1189 _ .oo81 ( ex P-) -0.0046 (theo.j _ . 055 (.Pdtj 

The quoted values are obtained for Q 2 > 
575 GeV , the lower Q 2 regions give consistent 
results. The experimental errors are dominated 
by the energy scale uncertainty and by the QCD 
model dependence of the detector correction. The 
theoretical uncertainty is given by the hadronisation 
correction (±0.0026) and the fi r dependence: 

+0.0053 r 2 . . n 2 ar .j +0.0025 r 2 . . F 2 

-0.0038 IOr Mr — V ancl -0.0007 1OT Mr ~ &T' 

The correlation between the fitted a s and the value 
assumed in the PDF seems to be - in particular at 
low Q 2 - stronger for jet rates than for jet cross 
sections. 



5. Summary 

a s (Mz) has been determined in NLO from jet 
observables in DIS at HERA (see Fig. |5j) . The result 
is consistent with the world average and has an error 
of about 0.006. This is quite a remarkable result, 
since it is only slightly less precise than the world 
average value having an error of 0.004 jy. 

In future, fitting technique using F 2 and jet data 
should be further developed to allow a simultaneous 
fit of a s and the PDF from HERA data alone. 
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Figure 5. Summary of recent ct s results from jets. 



References 

[1] S. Bethke, In IV th Int. Symp. on Rad. 

Corrections (1998), Barcelona (Spain) 

hep-ph/9812026. 
[2] T.Carli, hep-ph/9906541. 

[3] M. WOBISCH AND T. WENGLER, 

hep-ph/9907280. 
[4] S. D. Ellis and D. E. Soper, Phys. Rev. D 

48 (1993) 3160. 
[5] HI COLLAB., In XXIX Int. Conf. on HEP 

(1998), paper: 520, Vancouver (Canada), M. 

Wobisch, hep-ph/9811210. 
[6] HI COLLAB.,, In XXIX Int. Conf. on HEP 

(1998), paper: 533, Vancouver (Canada). 
[7] S. Catani and M.H. Seymour, Nucl. Phys. 

B 485 (1997) 291, Erratum-ibid. B 510 

(1997) 503. 

CTEQ collab., hep-ph/9903282. 

M. Wobisch, Thesis, in preparation (1999), 

University Aachen (Germany). 
[10] A.D. Martin et al., Eur. Phys. J C 4 

(1998). 
[11] J. Huston et al., Phys. Rev. D 58 (1998) 

114034. 
[12] V.D. Elvira, In XXXIII Rencontres de 

Moriond (1998), Les Arc (France). 
[13] CTEQ collab., Phys. Rev. D 55 

(1997)1280. 
[14] E. Tassi, Thesis, in preparation (1999), 

University Hamburg (Germany). 
[15] D. Graudenz, hep-ph/9710244v2. 
[16] S. Catani, Yu.L. Dokshitzer and B.R. 

Webber, Phys. Lett. 285 (1992) 291. 
[17] N. Tobien, Thesis, in preparation (1999), 

University Hamburg (Germany). 



